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Abstract 

Recent observations by the Fermi-hAT showed that there are delayed arrivals of 
GeV photons relative to the onset of MeV photons in some GRBs. In order to avoid 
a large optical depth, the minimal value of the Lorentz factor has been estimated to 
be higher than 1000 in some brightest bursts. In this paper, we present a detailed 
calculation of the time delay between the MeV and GeV photons in the framework of 
the magnetic-dominated jet model. We find that the time delay strongly depends on 
the saturated bulk Lorentz factor of the jet. Inspired by this fact, we use this model 
to calculate the Lorentz factors of the four brightest Fermi bursts. The results indicate 
that the Lorentz factors are much smaller than that obtained from the "single-zone" 
scenario. The short burst GRB 090510 has a minimal Lorentz factor 385, while the 
three long bursts GRB 080916c, GRB090902b and GRB 090926 have almost the same 
Lorentz factors, with an average value near 260. Another interesting result is that, for 
long bursts, GeV photons are emitted after the bulk Lorentz factor saturates. For the 
short GRB, however, MeV and GeV photons are emitted at the same phase, i.e., either 
in the expansion phase or in the coasting phase. 

Subject headings: gamma ray bursts: general - ISM: jets and outflows 



1. Introduction 



It is well known that photons with energy higher than mgC^ ~ 0.511 MeV in a local jet frame 
may annihilate into electron-positron pairs. The large optical depth of 77 annihilation, as well as 
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the Compton scattering, will restrain these photons from escaping the jet. However, the observed 
GRB spectra often peak in the MeV range and sometimes extend to the GeV range. This is the 
so-called "Compactness problem" . Since the optical depth is proportional to the inverse of the bulk 
Lorentz factor (F) of the jet, the "c ompactness problem" ca n be solved if we assume that the GRB 

The req uirement of the thin optical 

The 



outflow is moving with a large F (jReesI 119661 : 



Piran 



1999|) 



depth for the observed high energy photons sets a lower limit on T (jLithwick Sz Sari 



2001 



Burst and Transient Source Experiment (BATSE) data i ndicate that the Band-like spectra of most 
bursts have a thermal component at the prompt phase (|Rvddl2004l ). and th e measurement of the 
temperature allows us to constrain both T and the initial size rg of the flow (jPe'er et al.l 120071 ) . 



The investigation of GRBs has entered into a new epoch since the launch of the Fermi satel- 
lite in June 2008. The Fermi LAT instrument has observed several GRBs with photon energy 



as high as tens GeV. Within the frai nework of the simplified "single zone" mode 



tons set a very large lower limit on F (ILithwick &: Sari 



2004 



Granot et al. 



2008: 



20 id ) . For example, lAbdo et al 



Abdo et al 



2001 



Soderberg et al 



Ackermann et al 



20101 . 



2011 



2003 



20121 : 



the GeV pho- 



Razzaaue et al 



Ghisellini et al 



(2009 



outflow was Fr 



900. 



that Fr 



1200. 



2009al 

3) showed that the minimal Lorentz factor of GRB 080916c 



Ghisellini et al 



Ackermann et al.l (|2010l ) analyzed the spectra of GRB 090510 and showed 



(|2010l ) estimated the decelerating time of GeV emissions and ob- 
tained the Lorentz factor of GRB 090510 as large as 2000. However, an efficient physical mechanism 
to boost the outflow to such a large F is still unclear. 

An interesting featu re of the Fermi ob ser vations is that GeV photons often arrived seconds 



later than MeV photons (lAbdo et al 



2009a 



b c; Ackermann et al. 



2010 



2011 



). The explanation of 

this phenomenon s hould include both the intrinsic emission mechanism and the traveling process 
(jChang et al.ll2012l ). In some quantum gravity theories, photons can interact with the quantum 
fluctuation of the space-time, so high energy photons travel slower than low energy ones. Although 
this effect is very small, it can cause a detectable time difference after photons travel a cosmological 

Such Lore ntz invariance violation (LIV) ef- 



distance (Gambini et al. 



199S 



Elhs et al. 



20081 . 



2011 



fects lead to a natural time delay betwe en GeV and MeV phot ons (Schaefer 1999 ' Abdo 



et al 



2009c 



Boggs et al.ll2003l : iNemiroff et al.ll201ll ). As it was shown by lAbdo et al.l (|2009d ) and IChang et al 
(I2OI2I ). ■;he LIV effect is very small, we will neglect this effect in the following. 



Without considerin g the LIV effects, the delayed arrival of GeV pho to ns can also be explained 



by several GRB models (jPuran fc: Kumar 



2011 



Bosniak &: Kumar 



2011 



Duran fc Kumar (2011 



assumed that photons are emitted by electrons via synchrotron radiation, it takes more time for elec- 
trons to be accelerated to a large Lorentz factor in order to radiate GeV pho tons. iBosniak fc Kumar 
(l201in used the magnetic jet model, which was initially introduced by iDrenkhahnI (120021 ) and 
Drenkhahn &: SpruitI (|2002l ). to account for this phenomenon. According to this model, the optical 
depth is larger for high energy photons than that for low energy photons. GeV photons can only 
escape at a larger radius where the optical depth is below unity. 



In the magnetic reconnection model, the Band-type spectra can be produced from the photo- 
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sphere through th e magn etic dissipation (|Gianniosll2006l : iGiannios k. SpruitI 120071 : lGianniosll2008l ). 
Koers Gianniod (120071 ) first considered the neutron effects in this model, which was later used by 
Meszaros Reed (120111) to interpret the GeV time delay. I n the magnetic-domina ted but baryon- 
loaded model (jKoers Gianniod 120071 : iBeloborodovl I2OIII : iMeszaros k. Reed I2OIII ) , MeV photons 
can escape the plasma at the photosphere radius, which correspond to the prompt emission. How- 
ever, GeV photons are produced by the nuclear inelastic collisions between protons and neutrons 
at a larger r adius. In such two-zone scenario, the s trong constraint on the bulk Lorent z factor can 



be loo sened (jHascoet et al.ll201ll : IZhao et al.ll201ll : IZou et al.ll201lh . As pointed out bv lZhao et al 



(|201ll ). the optical depth depends not only on the energy but also on the emission angle. An 
average Lorentz factor rmin ~ 600 can be estir nated for GRB (3 80916c, GRB 090510 and GRB 
090902b in the two-zone model. In a similar way, IZou et al.l (j201ll ) assumed that the GeV photons 
were emitted at a larger radius than the MeV photons, and gave an analytical formula for Fmin by 
calculating the optical depth of a GeV photon going through the MeV photons shell. 



In this paper, we use th e magnetic-dominated jet model discussed by lKoers &: Gianniod (j2007l ) 
and IMeszaros &: Reed (j201ll ) to constrain the bulk Lorentz factor of GRB outflows. We show that 
the Lorentz factor of the short burst GRB 090510 can be as small as 385, while that of the three 
long bursts converge to about 260. The rest of the paper is organized as follows: In section [21 
we illustrate the magnetic-dominated jet model and the producing mechanism of MeV and GeV 
photons briefly. In section [3l we use the delayed arrival of GeV photons in four Fermi bursts to 
calculate the bulk Lorentz factor of GRB outflow. In section HI we discuss the validity of this 
model. Finally, conclusions are given in section [5l 



2. The magnetic-dominated jet model 



The hydrodynamics of the GRB outflow depends strongly on its geometry structure. In the 
magnetic-dominated jet model, the Lorentz factor of the outflow i ncreases with radiu s roughly as 
dPrenkhahnlbood : IPrenkhahn fc SDruitl[20o3 : iMetzger et aPboiol : iGranot et allboill i^ 



Tir) ~ 



V r > Tsat, 



(1) 



^The bulk Lorentz factor of the flow in the magn etic reconnection model originally took the form T{r) ~ 
??(r/rsat)^''^ for r < raat and r(r) « for r > rsat ijPrenkhahnl I2OO2I : IPrenkhahn fc SpruitI bood ). where they 
considered t hat the flow starts with t he Alfven speed at the initial radius ro. A compact form r(r) = (r/ro)^^'' 
was taken bv lKoers fc Gianniod (|2007l ). where ro is a length scale defined by specific com bination of the para meters. 



However, a Poynting jet can also be accelerated efficiently without reconnection process (iGranot et al 



20 If"), where 



r also takes the form F ~ {r /ro)^^'' , but ro denotes the width of the magnetic shell. Bosniak fc Kum ar (20f I) 
assumed that the format in Eq.([T| is valid at least in the interval of the Thomson- and pair- production-photosphere 
radii, and ro is roughly the same order of the radius where the jet is launched. Since it was unphysical for the jet to 
be accelerated to a large speed instantaneously, the Lorentz factor at the base ro was taken to be of order unity. In 
the present work, we take the idea of Bosnjak and Kumar, and write F as in Eq.(IT}. 
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where rg is assumed to be the base of the outflow, and rgat is the saturation radius, r/ denotes the 
ratio of the magnetic energy density to the baryon rest mass energy density at tq initiaUy. 

The injected baryons include both protons and neutrons. Initially, the neutron-proton jet 
accelerates as a single fluid where neutrons and protons have elastic collisions. When the n — p 
collision time scale is longer than the expansion time scale, the neutron component will coast with 
a terminal bulk Lorentz factor r„ at a characteristic radius, while the proton component is still 
accelerated. Thus, the neutron comp onent is embedded in a faster proton flow, and the jet becomes 
a compound flow (jBeloborodovlboilh . 



The cross section of n—p collision is a^uc ~ o"7r(c/urei), where 0",^ « 3 x 10~^^ cm^, and Vre\ is the 
relative speed of p to n. When frei c, the collision is inelastic. This occurs when the comoving 
expansion time tg^p ~ r/2cr becomes shorter than the comoving collision time ?a l/UpaT^c. 
Here n'p = Lx/Anr'^mpC^'qT is the comoving proton number density, L is the isotropic equivalent 
luminosity, and x = np/{np + rin) is the proton fraction of the baryon density. This gives the 
characteristic radius r^/ro = ij^x/l-qV"^, where ij^^ = {La.,^/ ^■nmpC^ro)^^^ — 1-32 x lO^L^''^"' 
Here we have adopted the Q = Qn x 10"" convention. Making use of Eq. ([1]), one obtains 



^54 ^0,7 



^'o \r]^x/2r]^ r > r^at- 

The pion production by the inelastic collisions is inevitable. A certain fraction of energy is 
carried away by neutrinos, which is an important prediction of the baryon loaded jet model. The 
vr'^ decay gives primary injected GeV photons. However, these photons undergo cascades and 
ca n not escape t he op aque jet. Interactions in the plasma are complex, more details can be found 



m 



Beloborodovl (|201lh . 



Suppose the final components in the jet contain photons with a Band-like spectrum, and the 
peak energy is around MeV. These photons start to be emitted when tt = n'pCTTr /2T ^ 1, which 
gives the Thomson photosphere radius, i.e., Tph/ro = T]j,/2r]T'^, where cjt ~ 6.65 x 10^^^ cm^ is the 
Thomson cross-section, and tjt = [Lax /4:'KmpC^roY^^ — 2.22 x IO^Lj^^Tq y''^. Using Eq. ([T|) for F, 
one obtains 

rph ^ f rj^{riT/2r]f/^ r < r^at, 
ro I r]^/2r]^ r > r^at- 

The simulation of magneto hydrodynamics (MHD) shows that the jet can form a conical struc- 



ture (jTchekhovskov et al.ll2008l ). After the jet exits the stellar envelope, the inner jet runs faster 
than the outer sheath. Thus, the Lorentz factor tapers off towards the edges. In such structure, 
the neutrons from the outer sheath can drift into the inner core. The relative radial Lorentz factor 
ratio between neutrons and baryons is larger than 1, which ensures that the collisions are inelastic 



(jMeszaros fc Reesll201ll ). Suppose the jet has an open angle 6, the transverse pion optical depth 
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can be expressed as r^r^x ^ n'paT^rOT = r]^{rQ/r){x6/r]). The jet becomes transversely optically thin 
(ttt^x = 1) at Ttt^x, which is defined as 



^x9 
Vn — - 



(4) 



The dynamical evolution of the jet depends on rj. If rj is large, for instance rj = 600r/600) the 
saturation radius rsat may be larger than r-j^, and r^,_L. Making use of F = (r/ro)^/'^, one obtains 
the following characteristic radii: 



~ 4.04 X IQi^Lgf rg/^^xo/g^r/goo/^ cm, 
rph ^ 3.98 X IQi^L^f r^Air/eoo''' cm, 
r-K,± ^ 4.41 X 10^^L54r/gQoXo.56'-2 cm, 
rsat =i 2.16 X lO^^rojr/loo cm. 



(5) 



On the other hand, if r] is small, rgat may become smaller than r7r,rph and r^r^^. In this case, one 
obtains the corresponding radii: 



(6) 



~ 1.32 X W^^L54Xo,5rilQf^ cm, 
rph ~ 5.98 X 10^^L54??foo cm, 
rn,± - 2.64 X 10^^L54r]~QQXo.59-2 cm, 
rsat ~ 1.00 X lO^^ro^rijlQQ cm. 

Here r] is in unit of 100, i.e., tj = 100??ioo- In both cases above, one has the order r^^ < rph < r.,^^±. 

The GeV photons are assumed to be produced at r^ _|_ by transverse nuclear collisions. However, 
the large optical depth prevents them to escape immediately. The spectra of produced photons 
depend on many parameters. Suppose that the Band spectrum peaks at -Ep ~ 1 M eV, then the 



optic al depth seen by a photon with energy E at radius r is approximately given by (jBeloborodov 



10^ 



40-/3-1 



r^2^L54 



E \ 



-13-1 



10 



GeV/ 



%00' 



(7) 



where (3 ~ —2.5 is the spectrum index above the peak energy Ep, and ri2 = r/10 cm. Setting 
r^^ = 1, we get the 7 — 7 transparency radius. 



' 77 



[E) K 2.50 X 10^^£;^/2r/go^oL54 cm. 



(8) 



Or equivalently, r^^{E) fa 1.94 x 10^'^ E'^^'^-q^^f^L^i cm for the rj ~ 100 case. Here, E is the observed 
photon energy in unit of GeV. Thus, at r7r,±, multi-GeV photons will be copiously produced by the 
transverse indrift neutrons colliding with jet core baryons. In the radius range r^ _|_ < r < r^^, these 
photons will annihilate into electron-positron pairs due to the large optical depth. Only beyond 
the radius r^^, the produced GeV photons can escape without obstructions. 
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The time delay for a photon with energy E relative to the onset time of MeV photons, equates 
to the time it takes for the jet to propagate from to r^^{E), 



At 



(1 + ^) / 



2r2( 



(9) 



where z is the redshift. The explicit formula for At depends on the order of rph, Tgat and r^^. From 
Eqs. (I5|), Q and ([8]), one can obtain rph < r^^ for the generic E > IQ GeV and r/ > 100. Using 
Eq. ([T]) for r, one obtains the formats of At for three different cases: 

( 3(l+^)ro 



At 



"77(S)-rsat I 



''sat ^ ^ph ^ ^77) 



(10) 



+ 



[l + z) 



rph < rs-at < r. 



77- 



A phenomenological illustration of the magnetic-dominated jet model is depicted in Fig. [TJ 
A magnetic-dominated but baryon-loaded jet is launched from a progenitor at the initial radius 
ro. The bulk Lorentz factor of the jet evolves as Eq.([T|). MeV photons are produced at the 
radius by nuclear collisions, vro-decay, electron-positron annihilation, magnetic dissipation, and 
synchrotron radiation, etc.. But these photons can only escape at the photosphere radius rph, 
where the Thomson optical depth decreases to bel ow unity. Thus, multi-M eV photons are emitted 
at Tph and lead to the observed Band-type spectra (jVeres fc Meszarosll2012l ). The GeV photons are 
assumed to produce at Ttt^j. by transverse drift nuclear collisions, inverse Compton radiation, etc.. 
But these GeV photons are capable to escape only at a larger radius r-y^ due to the large optical 
depth at r^,_L. The time it takes for the jet to propagate from rph to r^^ naturally leads to the 
GeV time delay relative to the onset of MeV photons. 



3. Constraints on the Lorentz factors 

From Eqs. dS]), (l6|), dH) and ([TO]) . the time delay At strongly depends on the terminal bulk 
Lorentz factor rj. Inspired by this fact, we make use of the magnetic-dominated jet model discussed 
above to calculate rj for four Fermi bursts, GRB 080916c, GRB 090510, GRB 090902b and GRB 
090926, respectively. 

The observed parameters which are necessary in the calculation are listed in Table [H Note 
that GRB 090510 is a short burst, while the other three are long bursts. In Table [H -Bhigh was 
taken to be the energy of the most energetic photon in each GRB. One exception is that the second 
energetic photon with -Ehigh = 11-16 GeV in GRB 090902b was chosen, while the most energetic 
33.4 GeV photon arriving at 82 s was excluded. This is because the isolated photon is far apart 
from other GeV photons and it is quite possible that this individual event happened when the jet 
encountered the interstellar medium. 



Fig. 1. — A phenomenological description of the magnetic-dominated jet model, rgat can locate 
before rph, or beyond r^^, or at any place between them. 



Table 1. 



GRB 


Eios,54 


Tgo 


z 




Atobs 


080916c 


8.8 


66 


4.35 


13.22 


12.94 


090510 


0.11 


0.6 


0.90 


31.0 


0.20 


090902b 


3.7 


22 


1.82 


11.16 


9.5 


090926 


2.2 


13 


2.11 


19.6 


21.5 



Note. — The observed parameters of four Fermi- 
detected GRBs. -Eios,54 is the isotropic equivalent en- 
ergy in unit of 10^^ ergs. Tgo is 90% the GRB du- 
ration time in unit of second, z is the GRB redshift. 
-E'high is the highest energy of photons for each burst in 
unit of GeV. Atobs is the observed time delay between 
the highest energy photon relative to the onset of 100 
MeV photons. The data were taken from Ichang et al] 
(I2OI2I 1. 
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There are several parameters which are uncertain, such as the initial radius rg, the jet open 
angle 6 and the ratio of proton number density to that of baryons x. Long bursts usually have 
time variability 6t ~ 10 ms, thus the initial radius is taken to be rg ~ c6t ^ 10^ cm for long bursts. 
The value of ro for short bursts is usually assumed to be smaller than that of the long bursts, and 
we set ro ~ 10'' cm for short burst GRB 090510. A nominal value of jet open angle 6 is taken to 
be 0.01, and the proton fraction of the baryon density x is approximately 0.5, i.e., 6-2 ~ 1 and 



j;o.5 ~ 1 (iMeszaros &: Reeal201ll ). \ yhen x approaches zero, this model reduces to the magnetic jet 



model without the loaded baryons ([Bosnjak &: Kumaijl201ll ) 



As was mentioned in section [2l the explicit formula for At depends on the order of rph, Tgat and 
r^^, which is not known previously. Thus, a self-consistent calculation should be taken carefully. 

Case I: First, we consider the rgat > > case (see the first formula in Eq. (|10p ). The cal- 
culated saturation bulk Lorentz factor and characteristic radii are listed in Table [2j The character- 
istic radii of the short burst GRB 090510 follow the order rgat > ^^77 > '''ph; which is self-consistent. 
However, for three long bursts, the results indicate rgat < ^^77; which are in contradiction with the 
assumption. The short burst GRB 090510 has a Lo rentz factor about 7 20, which is much lower 



than the prediction of one-zone models. For instance, lAbdo et al.l (j2009bl ) presented that the bulk 
Lorentz factor of GRB 090510 was as large as 1200. Our results indicate that GeV photons in the 
short burst are emitted before the Lorentz factor of the jet saturates. 



Case II: Then, we consider the rgat < ''ph < '''77 case (see the second formula in Eq. (fTO|) ). 
The results are given in Table [3j For all of the four bursts, we have _|_ < rp^. Thus, the 
transverse nuclear collisions happen inside the photosphere, and GeV photons are converted to the 
cascades. The overlap of the producing regimes of GeV and MeV photons is possible, but GeV 
photons are attenuated until r^^. The spectrum of GRB 090510 is fitted we ll by the Band fun ction 



plus a power-law component which dominates in the band above 30 MeV (jZhao et al.l 120111 ). this 
can be explained well by the Magnetic-dominated jet model. If these arguments are true, the bulk 
Lorentz factor of GRB 090510 is further reduced to 385. In Table El one also notice that for GRB 
080916c, Tph < ?"sat) which is not self-consistent. The bulk Lorentz factor of GRB 090902b and 
GRB 090926 in this case are calculated to be 245 and 252, respectively. 



Case III: Finally, we consider the rph < rg^t < ^^77 case (see the third formula in Eq. (fT0l) ). 
The results are listed in Table HI The data of GRB 090510 are absent, because any value of rj can 
not fit Atobs = 0.2 s by the formula. The minimal value of At is 0.34 s locating at ryeoo ~ 1- The 
data of three long bursts fit well in this case. The bulk Lorentz factors are 270, 252 and 258 for 
GRB 080916c, GRB 090902b and GRB 090926, respectively. 

So far, we have obtained the self-consistent Lorentz factors and characteristic radii for all of the 
four bursts, and summarized the Lorentz factors in Table [5j From the calculation above, the two 
conditions r^ < rph and r^ _|_ < r.y^ always hold. For short burst GRB 090510, the self-consistent 
cases are case I and II, and the corresponding Lorentz factors are 720 and 385, respectively. For 
the three long bursts, case III is valid for all of them, but they are all excluded in the case I. The 
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GRB 




rph/cm 


r^^/ cm 


^sat/cm 


rjr/cm 


r^,_L/cm 


080916c 
090510 
090902b 
090926 


0.39 
1.20 
0.32 
0.26 


5.24 X 10^3 
1.29 X 10^3 
6.80 X 10^3 

7.73 X 10^3 


1.78 X 10^*^ 
3.18 X 10^^ 
4.67 X lO^*^ 
3.09 X 


1.28 X 10^^ 
3.73 X 10^^ 
7.08 X lO^'' 
3.80 X lO^'' 


5.33 X 10^2 
1.31 X 10^^ 
6.90 X 10^2 
7.84 X 10^2 


1.51 X lO^'' 
6.74 X 10^^ 
2.32 X lO^'' 
2.87 X lO^'' 


Note. — The calculated saturation bulk Lorentz factors and characteristic radii under the 
assumption that rgat > ''77 > '^ph- We choose roj = 1 for short burst GRB 090510 and tqj = 10 
for other three long bursts. The characteristic radii of three long bursts are not self-consistent. 








Table 3. 








GRB 


??100 


rph/cm 


r-y^/ cm 


rsat/cm 


rjr/cm 


r^^x/cm 


080916c 
090510 
090902b 
090926 


2.58 
3.86 
2.45 
2.52 


1.37 X 10^^ 
6.31 X 10^5 
1.48 X 10^^ 
1.62 X 10^5 


1.09 X 10^^ 
7.16 X 10^^ 
1.38 X 10^^ 
2.80 X 10^*^ 


1.72 X 10^5 
5.75 X 10^4 
1.47 X 10^5 
1.60 X 10^5 


5.02 X 10^2 
1.91 X 10^2 
5.96 X 10^2 
5.88 X 10^2 


1.36 X 10^^ 
1.25 X 10^^ 
1.81 X lO^'' 
1.77 X lO^'' 


Note. — The calculated saturation bulk Lorentz factors and characteristic radii under the 
assumption that rgat < '^ph < f-yy- The parameters are the same as in Table [2j The characteristic 
radii of the GRB 080916c are not self-consistent. 








Table 4. 








GRB 


??600 


rph/cm 


r^^/ cm 


?'sat/cm 


r^r/cm 


r^^x/cm 


080916c 
090902b 
090926 


0.45 
0.42 
0.43 


4.82 X 10^3 
5.77 X 10^3 
5.71 X 10^3 


8.68 X 10^^ 
1.20 X 10^^ 
2.50 X 10^6 


1.97 X 10^5 
1.60 X 10^^ 
1.72 X 10^^ 


4.89 X 10^2 
5.86 X 10^2 
5.80 X 10^2 


1.31 X 10^^ 
1.77 X 10^^ 
1.74 X 10^^ 



Note. — The calculated saturation bulk Lorentz factor and characteristic radii for three 
long bursts under the assumption that rph < rgat < ''77- The parameters are the same as in 
Table [H 
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allowed Lorentz factors for the three long bursts seem to converge to an average value about 260. 
The saturation bulk Lorentz factor of the short burst GRB 090510 is reduced sharply, but still 
higher than that of long bursts. Another interesting feature is that GeV photons are emitted after 
the bulk Lorentz factor saturates for long bursts, i.e. in the coasting phase. For the short burst, 
however, both MeV and GeV photons are emitted either in the expansion phase or in the coasting 
phase. Case I can not happen for the long bursts, while case III can not happen for the short 
bursts. 



Discussion 



Besides the time delay, anoth er importa,nt feature of GeV emissions is that they last much 



longe r than the sub-MeV photons (|Gao et al 



200S 



Kumar 



20091 : 



Ghirlanda 



20101 : 



Ghisellini et al 



20ld ). For instance, the duration time of the sub- MeV photons is 55 seconds i n GRB 080916c, 
while photons with energy > 100 MeV last about 1400 seconds (jKumaii l2009l ). The observed 
decline of flux can be explained by the synchrotron radiation in the external shock (ES), i.e., 
oc t{3/3+2)/4j^/3/2 ^ _2_4 fo^ QR^g 080916c). The data of the initial 55 seconds are able to 
explain the observed X-ray and optical flux of the afterglow one day later. Thus, the GeV emissions 
have an afterglow origin. 

The spectrum and the light curve of th e GRB 090510 were also explained by the synchrotron 
radiation in the ES model (jGhirlandall2010l ). iGhisellini et al.l (j2010l ) studied the light curves of 11 
GRBs detected by LAT, and concluded that LAT fluxes decay in a common way cx t~^'^ for 
the four brightest GRBs studied in this paper. The LAT fluxes can be interpreted as the fireball 
emission in the radioactive regime. The spectra of the GeV emissions in some bursts showed a 
different power law from the Band function. Thus, the spectra and the light curves present strong 
evidences that GeV e missions have different origin with the sub- MeV emissions. As hinted by 



Ghisellini et al. 



(|20inl ). one can divide the "total emission time" of sub-MeV and GeV emissions 
: one is the overlap regime where both the sub-MeV and GeV photons are present; 



into two part^: 

another is the regime where only LAT photons exist. The latter can be named as the early afterglow. 

Notice that our calculation about GeV time delay is valid in the overlap regime, the GeV 
emissions in the early afterglow are not discussed. Once the outflow collides with the environment 
medium, the forward shock can also occur in the magnetic-dominated jet model. Both electrons 
and protons can be accelerated by the shock and form a power law spectral distribution. The 
characteristic frequency of the synchrotron radiation follows v = Tj'^pqB /27rmf>^pC, the synchrotron 
radiation of protons can be ignored compared to that of electrons. The produced photons have 
the same power law spectrum with electrons. Since the optical depth is small at so large radius. 



^X-ray and optical radiation usually arrive one day later, i.e., in the afterglow phase. We do not include them in 
the "total emission time". 
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photons are emitted immediately. The radiative fireball leads to the long duration time and light 
curve of GeV photons. Therefore, synchrotron radiation in the ES can explain the GeV emissions 
of the early afterglow in the magnetic-dominated jet model. 

Now we consider the spectra and light curve in the overlap regime. In the neutron rich 
environment, the i nelastic collisions p roduce pions, which further decay into photons with minimal 
energy TOT MeV ( Fan &: PiranI |2008|) . In the meanti me, the produced neutrinos can escape with 



the observed energy ~ O.IF GeV (lBeloborodovll201lh. Howeve r, these high energy neutrinos are 



difficult to be detected on the earth. iKoers &: Gianniod (|2007l ) estimated that less than 1 GRB 
neutrino event can be detected every year for nominal GRB parameters {z = 1). Thus, one can 
not exclude the baryon loaded model by the neutrino argument. If the spectrum of the protons is 
in power law, the resulted photons will also follow the same distribution. However, these original 
photons will quickly convert to via 77 reaction unless they are produced at a large radius 
where optical depth is below unity. The subsequent processes including Coulomb and Compton 
interactions are complex. Finally a Band like spectru m can form and the radiation becomes the 
observed photons in prompt phase ( Beloborodovl boiH ) . 



Wi thout the neutron component, the magnetic dominated outflow can dissipate energy effi- 



ciently (jPrenkhahn 



2002 



Drenkhahn &: Spruitll2002l ) . A non-thermal spectrum can be produced by 



the magnetic re c onnection model, a.nd th is spectrum is close to the observed prompt GRB emission 



( Giannios 



2006 



Giannios &: SpruitI 120071 ) . This means a broken power law (Band- like) spectrum 



can be produced with or witho ut baryons. The lumi nosity of GeV emissions in the overlap regime 
rose as L oc in most GRBs ( Ghisellini et al. 201ol'l. One exception is GRB 080916c, where the 
luminosity rose as L (x t^, which was a puzzle (jKumarl 12009 ). Our conclusion that the long bursts 
have the order < r^at < ^^77 may help to understand this puzzle. In this order, the jet is still 
in the expansion phase after the prompt emission, which means that F increases with time. If the 
jet is magnetic-dominated, one has F ^ (r/rp)^ and r = 2actT'^. In this way, /i = 1/4 leads to 



F oc t^/^. Since L oc t^F^ ( Ghisellini et al.ll2010l ). one can explain the puzzle. The light c u rve o f 
the GeV emissions in the overlap regime strongly favors the magnetic jet model. Gao et al. (j2009l ) 
also found that the physical composition of the GRB 080916c is likely magnetic. 

The spectra of GeV photons do not evolve with tir ne, and have a flatte r component (the 
slope intermediate between a and /3 of the Band function) (|Ghisellini et al.ll2010l ). These evidences 
strongly indicate that the GeV photons have a different p roducing mechanism. The two component 
GRB spectra were discussed bv I Veres Meszaroa (j2012l ) recently, where a dissipative photosphere 
gives the prompt MeV emission, while GeV emission are produced by the inverse Compton scat- 
tering. The model studied in this work also belongs to the two component case, many possible 
spectra are able to account for different GRBs. Therefore, the magnetic dominated jet model can 
explain many phenomenon of GRBs, such as the GeV time delay, the light curves and the spectra, 
etc.. 
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Conclusion 



In this paper, we have studied the bulk Lorentz factor of GRB outflow within the framework 
of magnetic-dominated jet model. We found that the emission mechanisms of the short and long 
GRBs are different. The long bursts have a unified bulk Lorentz factor around 260 in both case 
II and III. However, the Lorentz factor of the short burst is 720 in case I a nd 385 in case II . 
These values are much smaller than that obtained from the "one-zone" scenario. IZhao et al.l (120 111 ) 
calculated the Lorentz factor of GRB 080916c, GRB 090510 and GRB 090902b, and showed that 
r] ~ 600 could be consistent with observations in the "two-zone" scenario. Their values were still 
2 — 3 times larger than our result for long bursts. They also proposed that the Lorentz factor could 
be even lowered in the "multi-zone" scenario. The magnetic-dominated jet model discussed here 
is a kind of "multi-zone" model, according to which photons with higher energy are emitted at 
a larger radius. The L orentz factors for long bursts we obtained here were well inside the limits 
(~ 200 - 400) given bv lZou et all (|201lh . 



According to the magnetic-dominated jet model, the Lorentz factor depends on the initial 
radius ro, which is an undetectable parameter. In case I, only short burst GRB 090510 is self- 
consistent. As indicated in Figure El %oo is asymptotic to 1.7 when vqj goes to infinity. This 
means that the maximal Lorentz factor of GRB 090510 is about 1000 in this model. A small ro 
leads to a small tj. For instance, if r^j = 0.1, t] becomes 360. In case II, r/goo does not depend on 
ro,7, because the formulas of r-y-y and rph are independent of roj (see Eq. ([6]) and dH)). In case III, 
rjGoo depends weakly on roj, since r^^ are 2 ~ 3 orders of magnitude larger than rph- Thus, rph 
can be ignored in the calculation. The strong correlation between At and r] has an advantage: a 
small variation of At will not lead to a big change of r]. 

The time delay of GeV photons relative to the MeV photons can be well explained in the 
magnetic-dominated jet model. The bulk Lorentz factors of both long and short GRBs are reduced 
significantly. For GRB 090510, the possible minimal Lorentz factor is 385. For three long bursts 
GRB 080916c, 090902b and 090926, the Lorentz factors converge to about 260. The Lorentz factor 
of short burst is still larger than that of the long bursts. One common feature of the long bursts is 
that GeV photons are emitted after the bulk Lorentz factor saturates. In contrast, GeV photons 
in short burst can be emitted either in the expansion phase or in the coasting phase, and the bulk 
Lorentz factor in the former case is about one time larger than that in the latter case. 

The fact that three long bursts have a common Lorentz factor may imply that the long bursts 
have the same origin. One preva lent idea is tha, t long GRBs are caused by the co llapse of a massive 
star (such as Wolf-Rayet star) dwooslevl \l99^ : IPaczvhskil ^9981 : IWooslev et "aD [2OO6). The short 
duration time and the large Lorentz factor of the short burst may imply a different kind of central 
engine mechanism. For instance, short GRBs can originate from the inerger of two compact objects 



(such as NS -NS binary system and NS-BH binary system) (iGoodman 
ZhanJbood ). 



1986 



Meszaros &: Rees 



1992 
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Table 5. 





080916c 


090510 


090902b 


090926 


^ph ^ '"77 ^ ^sat 


X 


720 


X 


X 


'"sat ^ ''ph ^ ^77 


X 


385 


245 


252 


'"ph ^ '"sat ^ ''77 


270 


X 


252 


258 



Note. — The saturation bulk Lorentz factors for the four 
Fermz-detected bursts in three different cases, "x" denotes the 
inconsistent case. 




Fig. 2. — The relation between the saturation bulk Lorentz factor 77 and the initial radius tq for the short 
burst GRB 090510. 
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